Towards a complete, continuous, Wigner function for an ensemble of spins or qubits 
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I. INTRODUCTION 



The drive to realise quantum technologies has led to 
great advances in experimental capacity. Techniques 
such as quantum state tomography are rapidly becoming 
just as much part of the experimentalists toolbox as the 
theorists. One of the most useful examples of such a to- 
mographic method is the phase space distribution of the 
Wigner function [1 . One key element in this function's 
success is that it can be clearly linked to classical statis- 
tical mechanics and solutions to the Liouville equation, 
rendering in a clear and graphical manner the differences 
between classical and quantum states of a system. In- 
deed the Wigner function has been used to great effect 
in demonstrating the realisation of non-classical states of 
trapped light [2 . The usual application of the Wigner 
function is to systems that are expressed in terms of con- 
tinuous variables. However, there has been some signif- 
icant progress for finite dimensional systems where spin 
analogues of Wigner functions have been used to demon- 
strate entanglement in atom-chip circuits [3j. The usual 
approach here is to construct a pseudo probability distri- 
bution that maps information about the state of the sys- 
tem onto the surface of a sphere [4-9 . Such approaches 
have been very successful. For example, in [4 the authors 
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provide a method for constructing Wigner-like functions 
on the sphere's surface producing intuitively useful pic- 
tures for Fock states and spin coherent states. This func- 
tion can also be used to clearly demonstrate the non- 
classical nature of distinct superpositions of spin coherent 
states (so called spin-Schrodinger cat states) [10 . 

The Wigner function for systems of continuous vari- 
ables is unique up to an arbitrary overall phase factor. 
Moreover, the Wigner function satisfies a number of other 
very important properties, for example its marginal dis- 
tributions obtained through integration over a coordi- 
nate (e.g. position) provide the probability density in 
the canonically conjugate variable (e.g. momentum) |TI1 - 
13 . While the previously mentioned spin Wigner-like 
functions are impressively good at providing intuitively 
meaningful pictures they do not satisfy all the properties 
of a true Wigner function. For example, distribution as- 
sociated to the singlet state for two spins in [4^^ is every- 
where zero. Indeed these methods tend to only be good 
at analysing one total angular momentum subspace at a 
time (usually the largest one with, for N spin-half parti- 
cles, the angular momentum quantum number / = N/2). 
Hence, these phase space representations are not com- 
plete. 

Such shortcomings are known and there have been var- 
ious attempts to improve on these functions. These im- 
proved Wigner functions are usually defined on a discrete 
lattices or in higher dimensional phase spaces and tend 
to have most or all of the desired properties of Wigner 
functions [14-19 . However, there exists a quid pro quo 
and realising the full spectrum of properties of a Wigner 
function is done at the expense of not providing phys- 
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ically intuitive pictures. While mathematically of real 
value, in our view, it is also very desirable to produce a 
function where there are obvious parallels between the 
quantum state and its graphical representation as is the 
case for Wigner's original function [1 . 

In this paper we partially address this problem by in- 
troducing a new, easy to compute, spin Wigner-like func- 
tion that overcomes some of the difficulties found in [4]-[8j 
(such as being able to represent the singlet state for two 
spins). While not realising a complete representation of 
the state (and thus not being a true Wigner function), 
our function has many desirable features: it produces 
physically intuitive graphical representation; it can re- 
produce to good approximation previous Wigner-like dis- 
tributions; and it is useful for a larger subspace of states 
than previous methods. 



II. STRATEGY 

Our strategy for addressing this problem comprises 
several stages. First note that for N spins S'^ has eigen- 
values /(/ + 1) given in terms of angular momentum 
quantum numbers /= — 1,..., such that / > 0. 

We also have the operator Sz with quantum numbers 
m = — — / + 1, . . . , / — 1, /. As 6*^ and Sz commute they 
share common eigenvectors that are indexed by these 
quantum numbers. It is a textbook result that for > 2 
these vectors do not specify a basis, as {S^^Sz} is not 
a complete set of commuting observables. We can, how- 
ever, form a basis by introducing sufficient other observ- 
ables with quantum number (s) k which commute with 
Sx^Sy^Sz, such that {|k, /,m)} is an orthonormal basis. 
However, when considering general properties of angular 
momentum it is common to use the shorthand {|/,m)} 
leaving k to be specified later. We will adopt this con- 
vention here. 

Now we note that we can identify with the two indices 
/ and m with indices for two uncoupled harmonic oscil- 
lator Fock states |ni) and |n2). The idea being that, as 
it is straightforward to define Wigner functions for two 
harmonic oscillators, we can map the state of our system 
onto this Wigner function. 

In order to generate the Wigner function we choose, for 
each /, a particular value k/ for the additional observables 
k. A map ft from the A-spin Hilbert space to the Hilbert 
space for coupled harmonic oscillators is then defined by 

f2|k,^,m) = (l"^^^^ = l^ + "^ ^-"^^ ^ = ^^ (1) 
[0 kT^k/. 

The Wigner function W(^i,Pi, ^2,^2) for a density ma- 
trix p in its m representation is then 

^ {n[ n'2 1 9.p9) I ni n2 > lynm; (^1 , Pi ) Wn^n'^ (^2 , P2 ) , 

,77-2 ,ni ,n2 

where Wnn'{(liP) are Moyal functions pQH22] . which will 
be discussed in more detail later on in the text. 



While it may initially appear problematic that this 
Wigner function is four-dimensional it will turn out that, 
for many examples, we can make use of symmetry to re- 
duce this Wigner function to a three dimensional distri- 
bution that will well suit our needs. 

We note that W^(<7i,Pi, ^2,^2) is defined in terms of 
the density operator. In principle we are able to replace 
this with any operator. Whether or not it is possible 
that, as with the Weyl- Wigner formulation, an equivalent 
description of quantum mechanics can be given entirely 
in terms of an extension to our Wigner function remains 
an open question. What is certain is that this can be 
done for systems comprising one or two spins. 



III. CONNECTION TO SPIN 
REPRESENTATION 

The standard basis for the Hilbert space is 

\s1S2 ■■■sn) , 5a = t or I . (2) 

For A = 1 (i.e. one spin-half particle), we define and 
S± = Si ^ \S2 as follows [25] . 



53lt) = |m 53ii) = -ii;) 
s. It) = \\) s. \i) = s+ 

In the general case we set 



Si. 



(3) 



Si\siS2 - ' - sn) = \Si{si)s2- - - sn) ^\siSi{s2) " - sn) + 
...^\siS2---Si{sN)). (4) 

These satisfy the 5u(2) commutation relations, 

[Si.Sj] = icijkSk- (5) 

An eigenstate of {S'^^Ss) with eigenvalues (/(/ + 1),^) 
will be denoted \lm) (but note that for / < A/2 there 
may be many such eigenstates, more properly denoted 
|k,/,m)). 

IV. IMPLEMENTATION 

As stated in section [H] our approach to constructing 
Wigner functions is to relate the system of spins to the 
two-dimensional quantum harmonic oscillator. Recall 
that a basis \n1n2) for the Hilbert space can be con- 
structed from the ground state |00) using ladder opera- 
tors: 



\n1n2) := 



1 



The ladder operators satisfy the relations. 



(6) 



(7) 
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We will construct some new operators Ji , J2 , J3 using 
the Pauli matrices cr^ as follows: 



Ji = - ( a| al)ai 



(8) 



This may be thought of as an operator analogue of the 
Hopf fibration and is called the Jordan- Schwinger map. 
The reason for using this map to define Ji is that they 
will then satisfy the same commutation relations as the 
Si-. 



(9) 



It follows from Eq. ([8| that the basis states \n1n2) are 
also eigenstates for J3 and J'^ = -\- J2 -\- J^: 



J3 \n1n2) 
\n1n2) 



ni - 712 



2 

ni + n2 



nin2) 
ni + 712 



(10) 

+ 1 ) |nin2).(ll) 



In other words, the state \n1n2) has angular momentum 
eigenvalues / = (ni + n2)/2 and m = (ni — 77,2) /2. 



Thus we can now define a homomorphism by identify- 
ing a spin state's / and m eigenvalues with ni and 722 ac- 
cording to the above prescription. This homomorphism 
maps the A/'-spin Hilbert space to the Hilbert space for 
the 2D harmonic oscillator and identifies the actions of 
Si and Ji, in particular sending eigenstates of S'^^Ss to 
eigenstates of J^, J3. 

The Wigner function for any state in the A/'-spin 
Hilbert space will by definition be the Wigner func- 
tion of the associated state in the harmonic oscillator 
Hilbert space. The Wigner functions obtained will be 4- 
dimensional. We will show how this might work in prac- 
tice in the next section by considering some examples. 

For later use we briefiy recall the Wigner functions for 
the harmonic oscillator [13 . The Wigner function for the 
ID harmonic oscillator is defined by the formula 



(12) 



In the particular case where p = \n'){n\ one finds 



with denoting the Laguerre polynomials [24]. These 
are the Moyal functions for the harmonic oscillator [20- 
[22] , The Wigner function for the 2D harmonic oscillator 
is defined similarly: 

(14) 

where g{yi,y2) = {qi - yi,q2 - y2\p\qi +yi,q2+y2)- In 
the particular case where p = \n[n'2) {nin2\, 

1712 ,n^n2 iQl,Pl,q2,P2) = Wn^„'^iqi,Pl)Wn2„'^iq2,P2)- 

(15) 

We note that in use was also made of two harmonic 
oscillators to map \j = (ni + n2)/2,m = (ni — n2)/2) to 
\ni = j -\-m)\n2 = j — m), however our focus is different 
as these authors then go on to relate their work to that 

of H mill [9]. 

V. EXAMPLES 

We now construct linear maps Q from A^-spin Hilbert 
spaces to the 2D harmonic oscillator Hilbert space, and 
use these to construct Wigner functions. We require that 
the maps Q satisfy two conditions: 



(i) ns, \^) = j,n \^) 

(ii) rt^rt is a projection. 

Constructing Q in this manner is equivalent to construct- 
ing ft by introducing additional observables k, as de- 
scribed in section |ll| condition (i) is equivalent to de- 
manding that the additional observables commute with 
Si] and condition (ii) is equivalent to assuming that the 
states |k, /,m) are normalised. 

In order to make the Wigner functions as useful as 
possible, we choose the maps Q so as to maximise the 
rank of Q^ft and therefore the number of states that will 
be represented. 



A. One spin 

The 1-spin Hilbert space is 2-dimensional. The eigen- 
values of 5*^, 5*3 are summarised in the following table: 
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S3 


eigenspace dimension 


3/4 
3/4 


1/2 
-1/2 


1 
1 
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We will use these as a guide for constructing a map Q to 
the 2D harmonic oscillator Hilbert space. 

First, we observe that there is a unique state |t) with 
eigenvalues 3/4,1/2. We will map this to the unique 
harmonic oscillator state with the same eigenvalues: 



n|t) = |io). 



(16) 



As required by condition (ii), this choice preserves norms, 
since | t) = (10|10) = 1. We note that we could have 
multiplied the right hand side by a phase without spoiling 
this property. The action of Q on the other state is fixed 
by condition (i) above, via ladder operators: 



j-n\ 



(17) 



This implies that 0.\D = |01). Hence, the map that we 
have constructed is uniquely determined by conditions (i) 
and (ii) above, up to an unimportant overall phase. 

Having constructed the map f^, we compute some 
Wigner functions. The Wigner function for the state |t) 
is 

^io,io(<?i,^2,Pi,P2) = Wii{qi,pi)Woo{q2,P2)' (18) 

The Wigner function for the state |^) is similar, but with 
the roles of qi^pi and q2^P2 reversed. 



The Wigner function (18) is invariant under rotations 
of the form 



/«0 




/ 


COS^ 


sin 








\ 




Pi 






— sin^ 


cos 










pi 














cos 6 


sin 6 






\P2) 




\ 








— sinO 


cos 6 




\P2j 



(19) 

It follows that this Wigner function can be written as a 
function of xi^X2^xs G M^, defined using the Hopf fibra- 
tion by 



qi - m Q2 - m ^ 



q2 + ip2 



(20) 



Observe that the Wigner function of Eq. ( 18 ) is a map- 
ping from to R. A new function W : R^ ^ R may be 
constructed, making use of the above coordinate trans- 
form, in such a way that Eq. (18) can be rewritten in 
the form 



W(xi,X2,X3) = V^io,io(<?i,^2,Pi,P2) = -^e ''Li(x3+r) 

(21) 

(where r = ^y xl -\- -\- and we have used W instead 
of W to distinguish the function with transformed argu- 
ments form the original Wigner function). Isosurfaces for 
this function have been plotted in figure [l] 

We have also found it useful to represent our Wigner 
functions as functions on the sphere. Given a Wigner- 
like function W, we can integrate out the radial degree of 




FIG. 1: Isosurfaces of W for the state |t)- This Wigner- 
like function is positive within the blue regions and negative 
within the red regions. 



freedom and define another Wigner-like function of the 
surface of a sphere 

Ws(^,(/)) = - / W{x,y,z)rdr (22) 
4 Jo 

TT 

= — y W(r sin 6> cos 0, r sin 6> sin ^, r cos 6>)r dr. 

Although this may remove some information it does allow 
for direct comparison with preceding work |4J-19J . For this 
reason we have included, throughout this paper, plots of 
W together with those of W§. Note that the integral here 
is with respect to rdr, and not r^dr as one might expect 
from the volume form of R^, for the following reason. If 
the original state is normalised so that (V^IV^) = 1 
then the Wigner function on the sphere has integral 1: 



Is 



Ws((9,0) sini9di9d0 = 1. 



(23) 



The representation as functions on the sphere allows us 
to compare with the Bloch sphere representation: see 
figure [2j 

More generally, it can be shown that the Wigner func- 
tion W{qi^pi^q2,P2) for any operator that commutes 
with S"^ can be written as a function W on R^, and 
hence that the Wigner function for any eigenstate of S"^ 
is a function of three variables (see appendix [A|). On the 
1-spin Hilbert space S'^ has only one eigenvalue, so all 
Wigner functions can also be written as functions W on 



B. 2 spins 

The eigenvalues and eigenspaces of the operators 5*3, S'^ 
are: 
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eigenvalues (0,0) is again identified with a harmonic os- 
cihator state with the same eigenvalues: 






FIG. 2: Representation of the 1-spin state + \i))/V2 as a 
function in (top), a function on the sphere (middle), and a 
point on the Bloch sphere (bottom). Blue indicates positive, 
red indicates negative, and green indicates 0. 





^3 


eigenspace dimension 


2 


1 


1 


2 





1 


2 


-1 


1 








1 



As above, we use these to construct a map Q to the 2D 
harmonic oscillator Hilbert space. First, we identify the 
state Itt) with eigenvalues (2,1) with the unique har- 
monic oscillator state with the same angular momentum 
eigenvalues: 



^^^(in)-iit)) = |oo) 



(25) 



This completes the construction of the map Q. We 
now consider to what extent the map Q is unique. The 
only choices that we made in specifying Q were choices 
of phase in equations (24), (25). So there is apparently a 



2-parameter family of operators Q satisfying conditions 
(i) and (ii). Now one of these parameters corresponds to 
a choice of overall phase, so is unimportant. The other 
parameter can be compensated by making a rotation of 
the form ([l9|, so is again unimportant. Thus the map Q 
is uniquely determined by conditions (i) and (ii). 

Wigner functions can now easily be written down. 
Here are a few examples of Wigner functions for oper- 
ators that commute with 5^. These are all invariant un- 
der the rotation ([l9|, so can be written as functions on 
R^. Some of these Wigner functions have been plotted 
in figure [3j 



state 


w 


in) 

^(in) + i;t)) 

\U) 


:^e"''L2(r + X3) 
^e~''Li{r + X3)Li{r - X3) 
■^e~''L2{r - xs) 


7i(in)-i;t)) 





It is worth comparing our construction with those of [H 
E] at this point. Their construction yields interesting 
Wigner functions for the triplet states (i.e. eigenstates of 
S'^ with eigenvalue 2) but not for the singlet state, which 
is not defined. In contrast, our construction yields non- 
zero Wigner functions for all states, both the triplets and 
the singlet. 

For some states the Wigner function W is not invariant 
under the rotation of Eq. (19). We have not been able 



to determine a way to reduce such Wigner functions to 
functions W of three variables or functions W§ on the 
sphere. Indeed, Wigner functions for operators that do 
not commute with S'^ are not invariant under the rotation 



of Eq. ( 19 ) and we do not know if they can be reduced to 
sensible analogs of W or Wg. For example, the Wigner 
function for the operator |tt) ((til - (it|)/v^ is 



V2 



-Ql-Pl-(l2-P2 



(qi-ipif. (26) 



!^ in) = 1 20) 



(24) 



C. Three spins 



The action of Q on the other eigenstates of S'^ with eigen- 
value 2 is fixed by equation (17) above. The state with 



The eigenvalues and eigenspaces of the operators S^^ Ss 
are: 






FIG. 3: Isosurfaces of Wigner functions for the triplet states Itt)? (Iti) + lit)) /v^ and |ti), and the singlet state 
(Iti) ~ lit)) /v^- The Wigner functions are positive within the blue regions and negative within the red regions. 
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S3 


eigenspace dimension 


15/4 


3/2 


1 


15/4 


1/2 


1 


15/4 


-1/2 


1 


15/4 


-3/2 


1 


3/4 


1/2 


2 


3/4 


-1/2 


2 



Note that, unlike in the 1- and 2-spin Hilbert spaces, 
some of the eigenspace dimensions are greater than 1. 
This means that our method will be able to construct 
Wigner functions only for a subspace of the 3-spin Hilbert 
space. 

The construction of Wigner functions proceeds as 
above. First we fix the action of Q on states which are 
anihilated by 5+, guided by the eigenvalues of S'^ and 

^3: 



„27ri/3 



int) 



^(i;tt) + e-'"/'int) + e-^^^/^im)) = c|io). 



n itn) 

4.i/3 1^^^^^ 
47ri/3 I 



A 1 30) 

5 1 10) 



Here A, 5, C are complex parameters. The action of Q 
on the rest of the Hilbert space is fixed by equation (17). 

The choices made in the above three equations guar- 
antee that Q satisfies condition (i). Condition (ii) is sat- 
isfied only if the parameters A^B^C satisfy 



\A\ 



\B\ 



(27) 



These constraints reduce the number of parameters in 
the choice of Vt by two, leaving four parameters. One of 



these parameters corresponds to a choice of overall phase, 
and another corresponds to rotations of the form (19), so 
the total number of effective parameters is two. 

As in previous examples, the map Vt can be used to 
construct Wigner functions; we omit the details. 



D. N spins 

As the number of spins increases beyond 3, the number 
of parameters in the choice of ^ grows. However, the 
action of Vt on the "outer shell" (the eigenspace of S'^ 
with maximal eigenvalue) is always fixed, up to a phase. 
This is because the eigenspaces of 5*3 in the outer shell 
always have dimension 1. The action of on the outer 
shell is specified up to a phase by the equation 



1^ It . . . t) = \m) 



(28) 



and equation ( |17[ ). Here is a sample of Wigner func- 
tions obtained from states in the outer shell (where 
|(g) t) = It • • • t) and |(g) t) = It ... t) are spin coherent 
states [23 J: 



Operator 


W 


|(g) t) ((^ tl 
|(8) t) + ((^ tl 


i^e--L^(r + X3) 
^ J e ""LNir X3) 

^[(xi+iX2)^ + (xi-iX2)^] 



The results we present in Fig. |4] for spin Fock states 
bear striking resemblance to those found in [H [5] . In [10] 
the Wigner function of [H ^ was used to examine dis- 
tinct superpositions of spin coherent states (so called spin 
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Schrodinger cats). In Fig. [s] we present W and W§ for two 
example spin coherent states, the coherent superposition 
(cat) and their incoherent statistical mixture. In each 
case we see intuitively reasonable results that are in line 
with both our expectations and the results given by f4|l5l. 
Consider W and W§ for the cat and the statistical mix- 
ture shown in Fig. [5] Notice that the functions for the co- 
herent superposition of distinct spin coherent states (cat) 
contain oscillating fringes whilst those of the statistical 
mixture do not. As with Wigner's original function, the 
oscillating fringes indicate non-classical characteristics a 
quantum superposition of states. 

We note that, unlike the Wigner function presented in 
[HIS, our definition also allows for representation of some 
(but not all) states below the outer shell (for example, 
the singlet state for two spins as already discussed). A 
more detailed investigation is beyond the scope of this 
paper. We also note that the Wigner-like function W for 
spin coherent states has more than one local maximum. 
As such comparison to the Wigner function for coherent 
states of the harmonic oscillator must be done with care. 
However, for large A^, as with the function presented in [H 
[5], W§ becomes increasingly like the expected Gaussian 
distribution. 



E. Squeezed states 

As a final example application we consider the spin 
analogue of squeezed coherent states. Recall that 
squeezed states for the harmonic oscillator can be defined 
to be states of the form ex.p{f3{a^)'^ — ^o?) with \a) 
a coherent state and /3 G C a parameter. Analogously, 
we define squeezed states in the system of N spins to be 
states of the form 



exp(/?(5+)2-^(5_)2)|a;iV,0, 



(29) 



with \a]N^l) an spin coherent state [23 . Plots of the 
associated functions W and Wg are shown in figure 
[e] When is smah the Wi gner-like functions of the 
squeezed spin states resemble the Wigner functions of 
squeezed harmonic oscillator states, but for large val- 
ues of 1/3 1 the Wigner-like functions spread out and meet 
themselves on the opposite side of the sphere, leading to 
interesting interference effects. We note that these re- 
sults provide a good intuitive and accurate picture of the 
quantum state of the system. 



VI. CONCLUSIONS 

We have presented a new quasiprobability distribution 
function for ensembles of spins/qubits that has many 
properties in common with the Wigner function for sys- 
tems of continuous variables. Our function takes the 
form of a continuous distribution in which can, in 
many cases, be reduced via symmetry arguments to 'E? 



and the sphere. It enables clear graphical representa- 
tion of a wide variety of states, including Fock states, 
spin-coherent states, squeezed states, superpositions and 
statistical mixtures. Our definition respects the group of 
rotations generated by 81^82^ 83 just as Wigner's original 
function respects the action of the position and momen- 
tum operators. In addition, this property leads to our 
function exhibiting characteristics with clear parallels to 
other representations such as the Bloch sphere. More- 
over, even for relatively large it is easy to compute 
both W and Ws using standard mathematical software 
(to aid calculation of W§ analytical expressions are de- 
rived in appendix [B]) . 

Our definition allows us to go beyond previous at- 
tempts at constructing a continuous Wigner function for 
ensembles of spins [IHS!- Specifically we include states 
with different eigenvalues of 8'^. For example, with 2 
spins our construction yields a non-trivial Wigner func- 
tion for the singlet state. Nevertheless, not all states 
can be represented and for three or more spins/qubits 
information is still lost in passing from the Hilbert space 
to Wigner functions. The possibility of the existence of 
a continuous distribution that simultaneously is graphi- 
cally intuitive and forms a complete representation of the 
density matrix remains an open question. Nevertheless, 
we believe that our quasiprobability distribution provides 
a useful new tool for the analysis and visualisation of 
quantum states. 
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Appendix A: Rotational invariance of the Wigner 
function 



Our construction of the three-dimensional Wigner-like 
function W made essential use of the invariance of the 
4-dimensional Wigner function W associated to an eigen- 
state of 8'^ under rotations specified in Eq. (19). A short 
proof of this property follows. 



Let H'^^^ = ^(^^ +p^) denote the hamiltonian for the 
1-dimensional harmonic oscillator and let a denote any 
operator on the harmonic oscillator Hilbert space. Then 
the Wigner function for [ii^, a] is related to that for a as 
follows: 



(Al) 
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FIG. 5: Wigner functions for 5 spins. From left to right: the "up" state Ittttt)? the "down" state the "cat" (or GHZ) 

state Ittttt) + liiiii) and the statistical sum of the up and down states. The label G indicates a hidden Gaussian-like peak. 
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Suppose now that is an eigenstate of 5*^, or, equiva- 
lently, that the density matrix p = commutes with 

S'^ . Then VtpQ) commutes with . Now the hamiltonian 
H = ^1 + ^2 + i^i + 2-dimensional harmonic os- 

cihator is related to via the equation = — l)/4:. 
Since H has only positive eigenvalues, it follows that 
QpQ^ commutes with H. Then the 2-dimensional analog 
of Eq. ( |A1[ ) implies that 




FIG. 6: Wigner-like functions for squeezed states 
exp(/3(5'+)^ - ^{S-)^)\-tnn)- The left column shows 
^ = 0.1 and the right shows /3 = 0.2. The spherical function 
W§ has been plotted from two opposing viewpoints for 
greater clarity. 



Appendix B: Analytical expressions for 



Eq. (Al) may be proved by elementary operations: 



1 



{q-y\ [^q^a] \q ^ y) 



1 

^ J — oo 



d 1 
op TT 



{q-y\ \q^y)e' 



{Q-y\cr 



\q^y) 



i 



2 dqdy 



{{q-y\(T\q^y))e 



2\py 



d 1 

oq TT 



{{q-y\(j\q^y))e 



2\py 



When computing the spherical function W§ one en- 
counters integrals involving products of Laguerre poly- 
nomials and the exponential function. It will be shown 
below that these integrals can always be evaluated ana- 
lytically in terms of hypergeometric functions. Employ- 
ing these formulae results in a significant reduction in 
the time required to compute the Wigner function nu- 
merically. 



The spherical Wigner function associated to the oper- 
ator \lm) {lm'\ is 



4 7o 



l — m,l — m' 



{q2,P2)rdr 



{-ly ni+my.{i-m'y. C_sjn0ei*) 

47r Y {l+my.{l-m'y. \ '^^^'^^^ 



r 
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(Bl) 
(B2) 



where we have introduced 



POO 

irjic) := / e-Vi+"L^"((l + c)r)Lt{l - c)r) dr. (B3) 
Jo 

An integration by parts results in the equation 
Jo 

POO 

- {j + a) / e-V"L^_i((l + c)r)m{l - c)r) dr 



(i + a) J e-WL'^iil + c)r)L,"_i((l - c)r) dr. 

(B4) 



J 



Each of the integrals on the right of this equation eval- 
uates to a hypergeometric function [2j. The resulting 
expression simplifies on application of Gauss' relations 
for contiguous functions to 



X [(z + j + a + l)c^+^' + (j - i)c^+^-^] . (B5) 



While concise, this formula is numerically problematic 
since the hypergeometric function diverges at c = 0. 
Slight rearrangement yields the more practical formulae, 



(-1)'"' 1 + j + a; 1 + j - z; c^) [(i + j + a + l)c^-^ + (j - i^-'-^] i < j 



(i+g)! 



F{-j, 1 + i + a; 1 + i - j; c^) [(z + j + a + l^'^ + (j - i)c^-^-i] 



i > j. 



(B6) 
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